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Abstract 

We study the dynamics of molecular nanomagnets through a fully quantum mechanical model 
describing high-spin and high-anisotropy magnetic molecules subjected to a time-dependent 
magnetic field along the quantization axis, which continuously inverts the population of spin 
states. Crystals of molecular nanomagnets placed inside a resonant cavity interact with a 
quantized electromagnetic field. Relaxation of excited states takes place by means of spin-photon 
interaction, allowing stimulated emission of radiation and creating a maser effect. 
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I. INTRODUCTION 

Recent interest on high-spin high-anisotropy molecular nanomagnets has grown sig- 
nificantly due to the spectacular magnetic effects they exhibit, namely the pronounced mag- 
netic hysteresis and the quantization of the magnetic moment 1 . Potential applications of 
such magnets for current and future technologies include information storage, construction 
of nanomagnetic Maser-like devices, and quantum computation. The name Single-Molecule 
Magnets (SSM) has been coined to mean that individual molecules act as magnets. They 
can be prepared in long-lived excited quantum states by simply applying a magnetic field 2 , 
and exhibit a stepwise magnetic hysteresis in a time-dependent magnetic field. Quantum re- 
laxation of spin states at low temperatures is very sluggish, and has been intensively studied 
within the framework of the Landau-Zener-Stueckelberg(LZS) effect 3-5 . 

A detailed study of quantum spin transitions and the effect of the environment in an 
experiment with a rotating magnetic field, was done in Ref. 6. Their approach follows the 
lines of the LZS theory, taking as the starting point the following hamiltonian 6-8 : 



where D is the uniaxial anisotropy constant, g is the gyromagnetic factor, [i B is the Bohr 
magneton, B is the magnetic field and S is the molecular spin in units of ft. In the absence 
of an applied field B, the ground state is double degenerate, corresponding to the states 
with parallel (m = S) and anti-parallel (m = —S) projections of the magnetic moment 
along the quantization axis. Throughout this paper, we will consider only these two states 
(m = —S and m = +S) as relevant to the problem 6 . In this context, a striking effect is 
predicted 9 , where molecular nanomagnetic crystals exhibit a giant magnetic relaxation due 
to Dicke superradiance of electromagnetic waves 10 . In recent electron spin resonance (EPR) 
experiments, it has been observed a pronounceable resonant absorption of electromagnetic 
radiation by molecular nanomagnets 11-15 . In turn, they become a powerful source of coherent 
electromagnetic radiation when the wavelength of the emitted photons exceeds the linear 
size of crystals. When this condition is achieved, the molecules can coherently interact 
with the radiation emmitted, and the phase of the emitted photons may be considered the 
same throughout the sample 9 . Inside a resonant cavity, molecular magnets exhibit a strong 
dependence of the magnetization on the geometry of the cavity and this effect was observed 
experimentally, providing strong evidence for the coherent microwave radiation given off by 
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the crystals. These observations open the possibility of building nanomagnetic microwave 
lasers pumped by magnetic fields 16 . 

Our aim in this paper is to study the effect of stimulated radiation from molecular nano- 
magnetic crystals, by considering a fully quantum mechanical model, in which the elec- 
tromagnetic radiation enters as a quantized field. The essential ingredients in our model 
include: 

i) hamiltonian (1) describing individual molecules, with a time- varying magnetic field 
B z (t) along the quantization axis. This way, the low and high energy states will be 
continuously changing with time, creating the effect of spin-states population inversion. 
Note that this is different from most models proposed in the literature, where the 
magnetic field B z is constant, and the time-dependent field is applied in a transverse 
direction; 

ii) a quantized electromagnetic field inside a resonant cavity, which allows the relaxation 
of excited states by means of coherent photon emission, enhancing the radiation field. 

Finally, we will analize the semiclassical limit, in which the photon-field can be treated 
as a classical electromagnetic field. 

The above program will be developed in the present contribution. The content of this 
paper can be described as follows: in the next Section, we formulate the theoretical basis for 
analyzing spin dynamics, discussing the Hamiltonian. In Section III, we analyze the dynam- 
ics of a pure quantum mechanical state initially prepared in one of the double degenerate 
ground states of the system. We compute the correlation amplitude between the initial and 
the evolved state. Analytical approximated expressions are obtained to be compared with 
numerical results. In Section IV, we consider the photon field as a classical variable and 
obtain, in close analytical way, the conditions for a maser-like effect. Finally, in the last 
Section a few conclusions and remarks are added. 

II. HAMILTONIAN AND QUANTUM DYNAMICS 

We start considering the two lowest-level states (m = +S and m = —S) of spin 
Hamiltonian (1) and a second quantized term in photon variables describing the radiation 
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field in the cavity: 
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where the Pauli matrix 




yields the splitting of the levels in the presence of a time-dependent magnetic field B (t) 
applied along the quantization axis. The operator a^(a) is the photon creation (annihilation) 
operator. Inside a resonant cavity, only photons with a preselected frequency uj can be 
emitted or absorbed. The basis of Ho is given through the kets {\S, n >} and {| — S,n >}, 
being n the number of photons in a given state. We will consider that a general state can 
be written as follows: 

oo 

|tf >= {A S \S > +A_ S \ - S >} ® a n\ n > 

n=0 

The interaction between the molecular spin and the dipolar component of the electro- 
magnetic radiation inside the cavity is considered as a perturbation: 

/ 
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where the Pauli matrix 




flip the molecule from one level to the other (with the emission or absorption of a photon). 
One can see that our model consists of a two level system coupled to a quantized harmonic 
oscilator. Despite the simplicity of the above mentioned model, no exact analytic solution 
is yet known 17 . Now, we consider the quantum dynamics in the interaction picture, taking 
into account that the non-perturbed (diagonal) hamiltonian Ho is time dependent. In this 
case, operators evolve according to the equation: 

dO 

ih— = [O,H ] 
and for the interaction Hamiltonian we get 

Hj(t) = e iW ^H £ e- iW ^ , (4) 
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being 
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The temporal evolution of an initial state ket |\l/o > will be given by the unitary evolution 
operator Uj in the interaction picture. The latter can be written in the form of a Dyson 
series: 

Uj(t,0) = 1- % T /W iWr(tl) W/e- iW(tl) + 
n Jo 

+ ^ 2 j l dt^^Hie-^ J* dt 2 e m ^H ie - m ^ + ... , (6) 

with \ty(t) >= Uj(t, 0)|^o >• We can express the interaction hamiltonian Tij in the form 
below: 
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One must remember that H\ 2 and iJ 2 i are infinite-dimensional in the Fock space of photons 
We write Ui in the following form 

Ui = 1 + iF / dh 
Jo 
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The above expressions admit any time- varying magnetic field B (t). For sake of convenience, 



we will restrict our attention to the sinusoidal dependence: 

B (t) = B cos (Qt), (11) 

and in this case (8) and (9) reduce to 

H 12 (t) = a) exp [i (uot — r sin(fit))] + a exp [— % (cot + r sm(Qt))] , (12) 
H 21 (t) = a) exp [i (uot + r sin(fii))] + a exp [— % (uot — r sin(f2t))] , (13) 



where we have defined 
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The frequency V in Tii represents the coupling with the dipolar field, and is small when the 
wavelength of the radiation field is far longer than the molecular dimension. In this case, 
first order perturbation theory applies, leading to: 



Ui = 1 +iTMi + 0(T 2 ) 



and 



Mi = 



a* Fi + aF 2 
a^F* + aF* 

The dynamics of the system will be governed by the functions Fi(t) and F 2 (t): 
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being J2m+i(j') and J2m+2( r ) the Bessel functions of first kind, and 



Ai, m (t) = - we**" cos[(2m + 2) fit] - (2m + 2)VLe iuJt sin[(2m + 2)fii] (19) 
A 2 , m (f) = -iwe*** sin[(2m + + (2m + l)fie^* cos[(2m + - (2m + 1)Q. (20) 

/From the above formulae, a resonant behavior occurs when the photon frequency uj is an 
integer multiple of Q, say u = mil, and the dominant behavior of F\ and F 2 near that 
condition is given by: 

J m (r)[(-l) m i(l - e 2itjjt ) + 2tcj] 



Fx 



2u 



F 2 ~ (-1) F 1 = (-1) 



(21) 
(22) 



In the next Section we will analyze the dynamics of a pure quantum mechanical state initially 
prepared as an eigenstate of Hq. 
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III. DYNAMICS OF QUANTUM STATES 



To fix ideas, consider a quantum state prepared in the magnetic state \S > of the 
molecule, in the form: 

00 

|#o >= \S> E«nl»> ( 23 ) 

n=0 

where the mean number of photons at t = is given by 

00 

n =< era > = ^2 n \ a n\ 2 - 

n=0 

To first order perturbation theory, the evolved state \^(t) >= Ui\^ > reads 
\V(t) > « 1 x 

" " ! - 12 



X 



l + r 2 [(n + l)|F 2 | 2 + n |F 1 

^0 > +«r| - S >J2 a n (F*Vn+ l\n + 1 > +F*-v/n|n - 1 >) 



(24) 



n=0 

Observe that we have renormalized the evolved state. The correlation probability function 
defined as C(t) = | < ^ol^^) > | 2 > which is simply the probability of finding the system in 
the initial state at later times, yields 

C{t) = l + P[(no + l)|F 2 | 2 + no|Fi| 2 ] ' (25) 
and the above expression can be approximated at the resonant condition uo = mQ by the 
behavior of iq and F 2 given in (21) and (22): 

C(i) ^ i + HF,p ( 2„„ + ir (26) 

with 

|iq| 2 = |Jm(r , )|2 U sin 2 cut + Atcu sin 2ut + AcoH 2 ) . 
For long times ituj ^> 1), the leading behavior is quadratic in time 

|Fi| 2 ~ \J m (r)\H 2 , 
and we can approximate the above expression as follows: 

cm - , (27) 

with the correlation time r defined as: 

i = T| J m (r) I V2ti + 1. (28) 



The value 1/r can be interpreted as the rate of emission of photons, since the nanomagnet 
will relax to lower energy states by emitting photons. For large values of the constant r, 
maxima of the transition rate will be given by r = 2gSfiBBo/(hu/m) ~ (n, with (>la 
constant depending on the order of the Bessel function, yielding 



Next, we show some examples. To illustrate the evolution of a quantum state, we consider 
the initial mean number of photons no = 0, i.e. the vacuum for the cavity field, and the 
molecular spin in the S state (in this case the initial state is \^/q > = \S,0 >)• The resonant 
cavity is adjusted to the frequency / = uj/2-k = 10 GHz and the spin-photon interaction 
constant is taken as T = 2 GHz. In Fig. 1 to 3, we show the correlation probability C(t) 
as a function of time, for given values of r and m = u/Q. At a given ratio m = u/fl, 
the parameter r was set to give the maximum value of the Bessel function J m {r). The 
parameters for the first two figures satisfy the resonant condition (m an integer). One can 
see that the behavior for m even or odd is qualitatively the same. Except for the stepwise 
character of the correlation function, C(t) can be estimated by the expression (27), with the 
relaxation time given by (28). Note the irreversible character of C(t) at the resonance, due 
to the relaxation of the molecule by successive photon emissions. 

A quite different situation occurs for a non-integer value of the ratio u/Q, as shown in 
Fig. 3 for m = 9, 31. There are revivals of the initial state at latter times, in the form of an 
absorption-emission cycle not completely periodic. 

The above results were obtained within first order perturbation theory. Consequently, the 
long-time behavior should include corrections due to multiple photon processes. However, we 
believe that the essential physics is contained in the results discussed above. To support this 
view, in the next section we analyze the radiation field using the semiclassical approximation. 
We predict a masser effect at resonances. 

IV. SEMICLASSICAL THEORY 

The semiclassical theory is based on two fundamental assumptions: i) the cavity 
radiation field is described by a coherent state, which is the most nearly classical state (it 
minimizes the uncertainty relations); ii) the total density matrix is written as a product, 
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whose factors are related to spin and photon degrees of freedom. In other words, we assume 
that spin variables are uncorrelated with those of the field (Sargent) 18 . For a coherent state 
\z), we have the results 

a\z) = z\z), (29) 
(z\a* = (z\z*, (30) 
(n) = \z\ 2 . (31) 

So, the semiclassical approximation can be thought to be obtained by replacing the photon 
operators by complex numbers. Chosing the amplitude of the coherent state as real, we get 
a*, a — > yj (n). Then the equation of motion for the mean number of photons (n) will be 
given by: 

A A 

^ = -Tsm(ut)[p 12 e irsin ^ +p 21 e- irsin ^] 

with A = yj (n) The density matrix p of the spin system satisfies the following set of coupled 
equations: 

(U ' 4tT[p 21 e- irs[nm - Pl2 e irs[nm }Acos(ujt) (32) 



dt 

^ = -2ire~ irsin ^MA cos(cut) (33) 
<W = 2i Te irsin ^MA cos(^t), (34) 

where M = p n — p 22 . in general the solution must be accomplished by numerical methods. 
A simple solution is obtained if one assumes pu = p\\ = iN, being p± 2 — p 2 \ = 2iN, i.e. the 
off-diagonal elements depends on a single real parameter. It can be shown that the other 
possibilty, p\ 2 = p 2 \ = R does not lead to a maser effect. For the former case we have: 

dM 

— - = 8TNA cos(r sm(Qt)) cos(cjt), (35) 
dN 

— = -2TMA cos(r sin(Qi)) cos(wt), (36) 
dA 

— = 2TN sin(r sia(Qt)) emtut). (37) 
at 

An approximate solution near the values A ks 0, N and M for the photon field A is easily 
obtained yielding: 

A(t) =A + TN Re(F 1 - F 2 ) (38) 
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In the case of u = (2m + 1)0, the above solution for A grows linearly with time: 



A(t) « A + 2TNoJ 2m+ i(r) [t 



sin(2o;t) 



) 



(39) 



and consequently, the number of photons increases quadratically with time. In contrast, for 
uj = (2m + 2)fl the approximate solution is nearly constant, with A(t) ~ A , displaying 
ripples around this value. This even-odd symmetry breaking at resonance, is at variance 
with the pure quantum case treated in the previous section, where both instances presented 
similar behaviors at resonance, i.e. an increasing of the photon number due to relaxation 
of the molecular states, no matter if m was even or odd. In any case, for the semiclassical 
approximation, a maser effect is predicted for uj = (2m+l)f2. To illustrate this phenomenon, 
we show in Figures 4 to 7, the photon number as a function of time, for some temperatures. 
In Fig. 4 and 5, we take the limit ksT — > 0, and the initial photon number n — > 0(A — > 0). 
Other parameter values are M = and iVo = 0.5. In Fig. 4, we have chosen m = 9 and 
r = 10, while for Fig. 5, m = 10 and r = 12 (the corresponding value of r is always chosen 
to maximized the Bessel function). From Fig. 4, it is clear that the overall behavior of 
the photon number for m odd is closely related to the solution given by expression (39), 
and shown by the dashed line, corresponding to a quadratic increase of the number of 
photons with time. In contrast, when m is even, there is an oscillatory behavior, at least 
for short times, which reminds the quantum revivals. Figures 6 and 7 refer to examples at 
finite temperatures for m odd. For Fig. 6, T = 24 K, corresponding to the initial value 
no = kBT/(hu) = 49 for the number of photons, while Fig. 7 is at room temperature 
(T = 300K), corresponding to n = 625 photons. Note that the higher the initial number 
of photons, the faster the photon number increases with time. The effect is paramount in 
the latter case. Note also the different time scale. 

It practice, the photon number will increase until a saturation limit. The divergent 
behavior here obtained is due to the fact that losses are not taken into account. Such losses 
are provided by photons leaving the cavity and by excitation of other energy states. In 
fact, the real system which we are concerned here, is not as simple as a two-level object. In 
addition, the cavity field should include contributions from other photon modes, which we 
neglected from the beginning. 
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V. CONCLUSIONS 



We have developed the quantum dynamics of molecular nanomagnets, taking into account 
the quantized photon field of the cavity. In first approximation, the magnetic molecules were 
treated as two-level systems. The temporal evolution of quantum states were studied in de- 
tails in this paper, considering a time- varying magnetic field B z (t) = B cos(flt), applied 
along the quantization direction. The above field produces the necessary population inver- 
sion in a periodic way. 

The spin system couples with the dipolar component of the cavity field, allowing transi- 
tions between both molecular states with photon emission and absorption. The emission of 
photons is enhanced at the resonant condition, when the ratio between photon frequency uj 
and the applied field frequency Q is an integer number m = u/Q. At resonance, the energy 
is pumped from the applied magnetic field to the cavity radiation field, inducing a relaxation 
process depicted in Fig. 1 and 2. This irreversible process can be understood, if one con- 
siders that the probability of creating a photon increases with the number of photons. This 
striking phenomenon was illustrated in the previous section, via the semiclassical theory. 
In turn, the case of non-integer m produces a revival of the initial state, which means that 
emission and absorption of photons occur in almost periodic sequences, very similar to Rabi 
periods. 

Qualitative insights can be obtained via the semiclassical approximation, which is more 
suitable for practical purposes at room temperatures. Within this theory, the macroscopic 
state of the cavity field is described by the coherent states introduced by Glauber 19 . In this 
case, a maser effect is obtained at resonance, when m is an odd number. The maser effect 
is enhanced by temperature, when the initial photon number is increased. 
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Figure Captions 



Figure 1: Correlation function C(t) as a function of time, considering T = 2 x 10 9 s 1 , 
uj = 2n x 10 10 rad/s, Q = u/m, m = 8 and r = 9, 64. 

Figure 2: Correlation function C(t) with parameters T and u kept the same as for the 
previous figure, m = 11 and r = 11, 94. 

Figure 3: Correlation function C(t) for a non-integer ratio cu/fl. In this case we have 
chosen m — 9, 31 and r = 11, 45. 



Figure 4: Number of photons for the field A(t) = y (n) in the semiclassical approxi- 
mation using / = lu/2it = 10 GHz, T = 1 GHz, M = 0, N = 0.5, n = 0, m = 9 and r = 10. 



Figure 5: Number of photons for the field A(t) = y (n) in the semiclassical approximation 
using / = uj/2it = 10 GHz, T = 1 GHz, M = 0, N = 0.5, n = 0.25, m = 10 and r = 12. 



Figure 6: Number of photons for the field A(t) = y (n) in the semiclassical approximation 
with /, T, M = and N = 0.5 kept the same as the previous case, n = 49, m = 9 and 
r = 10. 



Figure 7: Number of photons for the field A(t) = y (n) in the semiclassical approximation 
using n = 625, m = 5 and r = 6.5. Other parameters kept the same as for the previous 
figures. 
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